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Abstract. Stellar catalogues in five passbands (UBVRI) over an area of approximately 0.3 deg , comprising about 
1200 objects, and in seven passbands (UBVRIJK) over approximately 0.1 deg 2 , comprising about 400 objects, 
in the direction of the Chandra Deep Field South are presented. The 90% completeness level of the number counts 
is reached at approximately U = 23.8, B = 24.0, V = 23.5, R = 23.0, / = 21.0, J = 20.5, K = 19.0. These multi- 
CIh band catalogues have been produced from publicly available, single passband catalogues. A scheme is presented to 

select point sources from these catalogues, by combining the SExtractor parameter CLASS_STAR from all available 
passbands. Probable QSOs and unresolved galaxies are identified by using the previously developed x 2 -technique 
^ | (Hatziminaoglou et al. 2002), that fits the overall spectral energy distributions to template spectra and determines 

the best fitting template. Approximately 15% of true galaxies are misclassified as stars by the % 2 method. The 
number of unresolved galaxies and QSOs identified by the x 2 -technique, allows us to estimate that the remaining 
level of contamination by such objects is at the level of 2.4% of the number of stars. The fraction of missing 
stars being incorrectly removed as QSOs or unresolved galaxies is estimated to be similar. The observed number 
counts, colour-magnitude diagrams, colour-colour diagrams and colour distributions are presented and, to judge 
the quality of the data, compared to simulations based on the predictions of a Galactic Model convolved with the 
estimated completeness functions and the error model used to describe the photometric errors of the data. By 
identifying outliers in specific colour-colour diagrams between data and model the level of contamination by QSOs 
is alternatively estimated to be < 6.3% in the seven passband and ^2.3% in the five passband catalogue. This, 
however, depends on the exact definition of an outlier and the accuracy of the representation of the colours by the 
simulations. The comparison of the colour-magnitude diagrams, colour-colour diagrams and colour distributions 
show, in general, a good agreement between observations and models at the level of less than 0.1 mag; the largest 
discrepancies being a colour shift in V — R and R — I of order of 0.2 mag possibly due to uncertainties in the 
bolometric corrections. Although no attempt is made to fit the model to the data, a comparison shows that the 
lognormal law for the initial mass function proposed by Chabrier (2001) describes the data better than the power 
law form in that paper. The resulting stellar catalogues and the objects identified as likely QSOs and unresolved 
galaxies with coordinates, observed magnitudes with errors and assigned spectral types by the x 2_ technique are 
presented and are publicly available. 

Key words, surveys - catalogs - methods: data analysis - quasars: general - white dwarfs - stars: low-mass, brown 
dwarfs 
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1. Introduction 

The primary goal of the ESO Imaging Survey project is 
to conduct public imaging surveys to provide data sets 
from which statistical samples of different types as well 
as of rare population objects can be drawn for follow- 
up observations with the VLT. A summary of the sur- 
veys conducted so far can be found in da Costa et al. 
(2001). Recently, deep optical/infrared observations of the 
Chandra Deep Field South (CDF-S) have been completed 
as part of the ongoing Deep Public Survey (DPS), and 
the data publicly released (Vandame et al. 2001; Arnouts 
et al. 2001). In these papers, fully calibrated images and 
single passband catalogues were presented for UBVRI ob- 
servations covering an area of 0.3 deg 2 and for UBVRIJK 
over a region of 0.1 deg 2 at the centre of the former area. 
Even though the angular coverage is still relatively small, 
the particular combination of depth and solid angle of the 
complete survey and the uniformity of the multi-band data 
which will result make it worthwhile to explore the type 
of objects likely to be sampled by the survey, to devise 
analysis techniques to identify objects of potential inter- 
est and to make a first evaluation of the samples that are 
likely to emerge when the full survey covering 3 deg 2 , is 
completed. 

The present work is in many ways an extension of that 
presented in Hatziminaoglou et al. (2002; hereafter H2002) 
which used a \ 2 -technique to classify point sources accord- 
ing to spectral type using their multi-band photometric 
properties, trying to identify special objects such as QSO, 
White Dwarf (WD) and low mass star candidates. 

Here, the question of how to choose an optimum stel- 
lar catalogue for statistical studies, either using a single 
passband or using multi-colour data, is discussed in more 
detail. Procedures are developed to set the contamination 
by extended sources below a specified value at the sin- 
gle passband level. To reach fainter magnitude limits the 
leverage of using multi-colour data is investigated. Finally, 
to assess to some extent the quality of the final multi- 
colour stellar catalogue produced, number counts, colour 
distribution, and colour-colour and colour-magnitude di- 
agrams are compared to those predicted by a Galactic 
model calibrated using independent data. 

In Section ^ the data are briefly reviewed as well as 
the method employed in the construction of the point 
source catalogue used in the present paper. Section || dis- 
cusses how point sources are selected in the single pass- 
band, and colour catalogues, and how potential unresolved 
galaxies and QSOs are found using a previously developed 
X 2 -technique. In Section [| the actual procedure for the 
creation of the single passband to the final stellar colour 
catalogue is presented. Section |E| and Section || introduce 
the Galactic Model and the simulation of appropriate data 
sets. Section presents the results. Finally, in Section || a 
brief summary is presented. 



2. The data 

The single passband catalogues for the CDF-S are taken 
from Vandame et al. (2001) and Arnouts et al. (2001), cut 
at a S/N — 2, slightly lower than the catalogues publicly 
released which were cut at a S/N = 3. 

This is done as a magnitude, even with a larger error 
bar, is preferred over an upper limit which would otherwise 
be assigned to the magnitude in a filter where an object 
is not detected. 

The £/-band catalogue used here (and in H2002) was 
extracted from an image produced by stacking all the 
available U/38 and C/350/60 images, to reach a fainter 
limiting magnitude. A more detailed discussion of this cat- 
alogue and of its photometric calibration will be present 
in Arnouts et al. (in preparation). The main conclusion 
is that a shift of +0.2 mag has to be applied to the U- 
magnitudesQ given by SExtractor running on the stacked 
U image using the photometric zero point of the U /38 
image, to get agreement between observed and predicted 
U — B colours for normal stars. 

To build a colour catalogue, the single passband cat- 
alogues are associated on position using, in the present 
paper, a fixed search radius of 0.8". With this search ra- 
dius the number of multiple associated sources is kept low 
(<,0.8% in all passbands, and typically 0.2%), while result- 
ing in many positive matches. A detailed discussion on the 
optimum search radius, and more complicated schemes 
using variable search radii will be presented elsewhere 
(Benoist et al. 2002). 

After the association, only objects that are in the area 
common to all catalogues and outside all the masks placed 
around saturated objects and bright stars are kept. This 
ensures that all objects could, at least in principle, have 
been detected in all passbands. The effective area is 0.263 
deg 2 for the five passband catalogue and 0.0927 deg 2 for 
the seven passband catalogue. If an object is not associ- 
ated in a passband, an 3<r upper limit is assigned for the 
magnitude in that passband. 

3. Selecting Point Sources 

3.1. Selecting Point Sources at the single passband 
level 

The SExtractor software package used at the moment to 
detect and extract sources uses a neural network to sep- 
arate point from extended sources. For each object it re- 
turns a parameter, CLASS_STAR, with a value between 1 (a 
perfect point like object) and 0. An example of the distri- 
bution of CLASS_STAR as a function of magnitude is shown 
in Fig. |l| for the i?-band. 

The neural network file used here is the standard one 
that is part of the SExtractor package and that has been 

1 All magnitudes in the present paper are in the natural sys- 
tem of the WFI filters at the ESO/2.2m telescope for UBVRI, 
and the SOFI filters at the ESO/NTT for JK, with Vega based 
zero points. 
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trained for seeing FWHM values from 0.025 to 5.5" and 
for images that have 1.5 < FWHM < 5 pixels, and is opti- 
mised for Moffat profiles with 2< (3 < 4. These conditions 
are fulfilled for the EIS optical and IR observations of the 
CDF-S. 

At the bright end of the distribution shown in 
Fig. ^ two sequences can easily be distinguished with 
class_star values below 0.1 and above 0.95, respectively. 
However, for fainter magnitudes the sequences spread out 
and merge. In previous EIS papers (e.g. Prandoni et al. 
1999) the notion of a morphological classification limit 
MAG_STAR_LIM was introduced, where objects brighter 
than this limit and with a CLASS_STAR value larger than 
a limit CLASS_STAR_LIM were considered point sources. In 
previous papers, these limits were chosen basically by eye 
from figures like Fig. [j]. 




20 25 
R— Magnitude 



Fig. 1. SExtractor parameter CLASSJ3TAR versus magni- 
tude for the i?-band, with the finally adopted morphologi- 
cal classification limits indicated by the solid line. Objects 
in the upper left box are considered point sources in a sin- 
gle passband. The dashed line indicates the much simpler 
criterion used previously which only took into account the 
CLASS_STAR value at the faintest MAG_STAR_LIM. 

Clearly, there is a need to determine CLASS_STAR_lim 
and MAG_STAR_LIM in a more systematic way, by statis- 
tically describing the relative contribution of the point- 
and extended-source population. The goal is to produce a 
magnitude-dependent classification as illustrated in Fig. ||. 

A technical description of the actual implementation 
of the concept used in the present paper is given in the 
Appendix, and Fig. [I] shows for the i?-band the dis- 
tribution over CLASSJ3TAR and the adopted values for 
CLASS_STAR_lim and mag_STAR_lim (solid line) in the 
present paper and for comparison the old way of using 
a single value (dashed line). It should be noted that this 
procedure allows the estimate of the smallest CLASS_STAR 
where stars can be distinguished over the large back- 
ground of galaxies. It does not imply, especially at the 
faint end, that stars with CLASS_STAR values lower than 
CLASS_STAR_LIM do not exist. Furthermore, it is recalled 
that the actual distribution of CLASS_STAR versus mag- 



nitude depends on the seeing but also on the Galactic 
coordinates of the field under consideration. 

It is important to emphasise that more work should be 
done to explore alternative ways of morphologically clas- 
sifying images. One alternative under consideration is to 
use multi-resolution methods to provide an independent 
pixel-based classification scheme (Vandame, private com- 
munication) . 

Apart from a statistical analysis of the distribution 
of CLASS_STAR versus magnitude a second criterion of se- 
lecting point sources is investigated, namely the fact that 
they are expected to be "round" . This is done by consid- 
ering the SExtractor parameters A and B which represent 
the profile fitting rms along the major and minor axis (see 
Arnouts et al. 2001 for more discussion). Visual inspection 
of some objects with large values of | ( A — B)/ (A + B) \ 
showed that these outliers are not stellar, but bad pixels, 
or objects close to the edge of a frame or a mask. The 
standard deviation of the distribution is computed and 
objects that deviate by more than a specified level (in 
the present paper 6a) are not considered point sources. 
Typically, 0.1% of the objects are removed in this way. 
A 3cr clipping, however was verified to be too stringent, 
as then many true stars would be removed. Most of them 
would be recovered by the colour Point/Extended (here- 
after, P/E) source classification scheme discussed below, 
but they would, incorrectly, be excluded from the single 
passband number counts. 

3.2. Selecting Point Sources in a colour catalogue 

The (P/E) source classification in a colour catalogue 
is a non trivial matter, as the MAG_STAR_LIM and 
CLASS_STAR_lim of all the passbands involved must be 
combined in some way. Here the following scheme is 
adopted. As the (P/E) source classification is most re- 
liable for bright objects (see Fig. |l|), the classification in 
the colour catalogue will be determined by the passband 
in which the object is brightest relative to the (faintest) 
MAG_STAR_LIM in that passband. If the object is classified 
as a point source in this passband, it is considered a point 
source in the colour catalogue. The fact that the seeing in 
the images is different is implicitly taken into account by 
the fact that the mag_star_lim in a case of worse seeing 
is brighter than in the case of good seeing. 

A completely different way of creating a colour cat- 
alogue is by digitally stacking the images in the various 
filters, and then running the source detection algorithm 
on it. The photometry is then obtained by going back 
to the individual filters and measure the flux around the 
now known positions. This method was applied to EIS 
data by da Costa (1998) using the algorithm of Szalay et 
al. (1999). From that work, SExtractor's standard neural 
network for (P/E) source classification is known not work 
properly on such a stacked image because of e.g. seeing 
variations among the different filters. Therefore this pro- 
cedure is not further pursued here. 
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3.3. Eliminating non-stellar objects 

The procedures outlined above can only separate point- 
like from extended objects. However, to produce a truly 
stellar catalogue one has to consider the possibility of 
contamination by QSOs and unresolved galaxies. This 
is not possible in a single passband catalogue. However, 
in a multi-colour catalogue the photometric data pro- 
vides information on the shape of the Spectral Energy 
Distribution (SED) which can, in principle, be used to 
address this issue. 

In fact, this was the subject of the paper of H2002 
where a colour point source catalogue was used to classify 
sources according to their spectral type. This was done by 
X 2 -fitting template spectra to the observed magnitudes 
in different bands. The spectral library in use consists 
of series of: model QSO, WD and brown dwarf spectra; 
three empirical cool WD spectra; a set of stellar tem- 
plates for O-M stars; and a set of model galaxy spectra. In 
H2002 emphasis was given in identifying QSOs, WD and 
low mass star candidates for follow-up spectroscopy. The 
X 2 -technique provides an independent view of the nature 
of a source, complementing the morphological classifica- 
tion. However, by using broadband photometric data to 
assign spectral types, misclassifications are possible (e.g. 
due to the fact that an object is an unresolved binary lead- 
ing to composite colours) and this is further discussed in 
Sect. Wl. 



Table 2. Number of objects in the five and seven passband 
point source catalogues 

UBVRI 



3. 4. Choosing parameters 

The morphological classification limits adopted in the 
present paper are listed in Tab. |l| and are chosen after 
extensive testing. In terms of the technical description in 
the Appendix the parameter £2 is set to 30 and the pa- 
rameter t\ to 80. 

The main guideline to chose these particular values for 
t\ and t 2 to set the mag_star_lim and class_star_lim 
values is the estimate of the remaining contamination by 
QSOs and unresolved galaxies, and related, the incor- 
rect removal of true stars by misclassification of the \ 2 ~ 
technique, which one wants to keep small (of order a few 
percent at most). These numbers are estimated indepen- 
dently in Sects. 17.2.2] and ffJ5. 



4. The procedure 

For the chosen combination of CLASS_STAR_lim and 
MAG_STAR_LIM, the single passband catalogues are pro- 
duced and then associated to form the colour catalogue 
with the (P/E) source classification applied based on 
all the available colours. This is done separately for the 
U BVRI and U BVRIJK catalogues. The resulting num- 
ber of point sources is given as the first entry in Table g. 
The second entry gives the split over the passband that 
actually decided that a source is considered a point source. 
In the five passband catalogue this is mainly R, I and B; 
in seven passbands mainly _R, J and B. 



Number of point sources 


1371 


Passband that decided (UBVRI) 


7/308/0/650/406 


Not classified by x 2 method 


18 


Classified by x 2 method 


1353 


as rank 1 


947 


as rank 2 


127 


as rank 3 


279 


Number of unresolved galaxies 


37 


Number of QSOs 


136 


Number of stars 


1198 


UBVRIJK 


Number of point sources 


486 


Passband that decided (UBVRIJK) 


2/110/0/225/37/110/2 


Not classified by \ 2 method 


5 


Classified by x 2 method 


481 


as rank 1 


160 


as rank 2 


117 


as rank 3 


204 


Number of unresolved galaxies 


16 


Number of QSOs 


61 


Number of stars 


409 



As a second step the x 2 template fitting of H2002 is ap- 
plied to the objects classified as point sources. This is done 
twice: once with and once without the galaxy template 
spectra included. The outcome of the template fitting is 
the spectral type of the best fitting template spectrum 
(with associated photometric redshift in the case of QSOs 
and galaxies), the (reduced) x 2 an d the "rank" (which 
is based on the x 2 value, rank 1 being the robust candi- 
dates, rank 3 poor candidates). Objects that are ranked 3 
when no galaxy templates are considered (typically hav- 
ing reduced x 2 larger than 4-7), and become ranked 1 
galaxies when galaxy templates are considered, that are 
fainter than 16th magnitude (in the passband where they 
are brightest) and with photometric redshifts larger than 
0.25 are considered to be (unresolved) galaxies and are re- 
moved from the list of point sources. The split over the as- 
signed ranks, the number of unresolved galaxies and QSO, 
and the number of objects considered true stars are given 
in Table |. 

The x 2 -technique requires an object to be detected in 
at least three passbands. Less than 1.5% of objects does 
not fulfil this criterion (see the third entry in Tab. |2|). They 
are kept as stars, as they are originally classified as point 
sources and there is no evidence to the contrary. Exclusion 
of this small number of objects does not affect any of the 
conclusions of this paper. 

From an observational point of view, the validity of 
the adopted approach to remove the most likely galaxy 
and QSO interlopers can be verified only by spectroscopic 
follow-up of a representative sample of objects. From a 
theoretical point of view it is planned in the future to 
simulate the physical size, redshift and colour space occu- 
pied by galaxies and QSOs as well, so it can be verified 
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Table 1. Morphological classification limits MAG_STAR_LIM (m x ) and CLASS_STAR_LIM (c x ) adopted in the single 
passband catalogues (see Fig. || for a visual representation). 



Passband 


Cl 


mi 


C2 


iri2 


C3 


m 3 


C4 


m,4 


U 






0.851 


21.57 


0.873 


22.27 


0.952 


22.45 


B 






0.848 


23.04 


0.933 


23.75 


0.961 


23.84 


V 


0.590 


21.15 


0.828 


22.15 


0.879 


22.68 


0.935 


22.69 


R 


0.863 


20.98 


0.873 


21.89 


0.907 


22.41 


0.915 


22.61 


I 


0.603 


19.30 


0.663 


20.25 


0.744 


20.77 


0.885 


21.09 


J 






0.644 


19.32 


0.929 


20.18 


0.923 


20.36 


K 










0.750 


17.71 


0.943 


18.46 



if a unresolved source within a given colour bin is more 
likely to be a star, galaxy or a QSO. 

The final number of objects considered to be stars is 
about 1200 over the 0.263 deg 2 for the five passband cat- 
alogue and about 400 over the 0.0927 deg 2 for the seven 
passband catalogue. For comparison, the number of ex- 
tended objects over these areas is about 74000 and 28000, 
respectively. 



5. The Galactic model 

The observed counts and colours will be compared to a 
newly developed Galactic model (Girardi et al. in prepa- 
ration). In brief, this model is based on a population syn- 
thesis approach, where stars are randomly generated from 
an extended library of evolutionary tracks according to 
user defined star formation rate (SFR), age-metallicity re- 
lation (AMR), and initial mass function (IMF), and then 
distributed in space according to the probabilities given 
by the specified geometry. The evolutionary tracks include 
brown dwarfs down to 0.01 Mq, and evolved phases with 
the complete white dwarf stage for stars up to 7 M Q . 
For the present paper the following ingredients have been 
adopted: 

— The IMF is taken to be the lognormal form proposed 
by Chabrier (2001). The effect of adopting the power 
law form proposed by him will be detailed in Sect. 7.1.4 

— The SFR is assumed to be constant over the last 11 
Gyr for the disk, and constant between 12 and 13 Gyr 
for the halo. 

— The AMR for the disk is taken from Rocha-Pinto et 
al. (2000). [Fe/H] values are converted into the metal 
content Z by means of a relation that allows for a- 
enhancement at decreasing [Fe/H], as suggested by 
Fuhrmann (1998) data. At any age, [Fe/H] is assumed 
to have a la dispersion of 0.2 dex. 

The metallicity of the halo stars is assumed be Z = 
0.0095, with a dispersion of 1.0 dex. This is based on an 
observed [Fe/H] value of — 1.6±1.0 (Henry & Worthey 
1999), allowing for an a-enhancement of 0.3 dex. 

— The disc component is described by a double- 
exponential in scale height and Galactocentric dis- 
tance. The model does not have separate components 
for thin or thick disk. Instead the scale height for disk 



stars is a function of age, and is parameterised as: 

H = z (l + t/t ) a (1) 

with t the age of the star. Rana & Basu (1992) for ex- 
ample find z — 95 pc, to = 0.5 Gyr and a = (2/3). 
However this does not fit very well the derived scale 
height of 'thick', 'old', 'intermediate' and 'young' disk 
components as derived by Ng et al. (1997). Their re- 
sults are described by zq = 95 pc, to = 4.4 Gyr and a 
= 1.66, which are adopted here. 

For the present paper the (relative) local densities 
of halo and disk stars, the oblateness of the halo and 
the Galactocentric disc scale length are derived from 
the "Deep Multicolor Survey" (DMS; Osmer et al. 1998 
and references therein) covering six fields and EIS data 
(Prandoni et al. 1999) for the South Galactic Pole. A 
more detailed discussion of the calibration of the Galactic 
model will be given in Girardi et al. (in preparation). Note 
that the Galactic model parameters derived below differ 
slightly from those used by H2002 to compare some of 
their results to. 

The halo oblateness (and local halo number density) 
is derived by fitting the number of halo stars, defined by 
(in the Johnson-Cousins system) < B — V < 0.7 in 
the range 20 < B < 22 and < V - I < 0.8 in the 
range 18 < / < 20, in these seven fields and is found to 
be q = 0.65 ± 0.05. This value is smaller than the value 
of 0.8 ± 0.05 quoted by Reid & Majewski (1993), but 
recently Robin et al. (2000) could not exclude a spheroid 
with a flattening as small as q = 0.6 and Chen et al. (2001) 
derived q = 0.55 ± 0.06. 

The disc radial scale length (and local disk number 
density) is derived by fitting the number of disk stars, 
defined by 1.3 < B - V < 2.0 in the range 20 < B < 22 
and 1.8 < V — I < 4.0 in the range 18 < J < 20, and is 
found to be TJr = 2800±250 pc. This is in agreement with 
the lower limit of 2.5 kpc (Bahcall & Soneira 1984) and 
the recent work of Zheng et al. (2001) on M-dwarfs who 
derive H R = 2750 ± 160 pc and Ojha (2001) who derive 
H R = 2800 ± 300 pc for the thin disc. 

The model with these parameters is applied to the 
CDF-S field. The Galactic bulge has no contribution in 
this direction. The model is calculated for an area of 3 
deg 2 , which is about 10 times larger than the actual area 
covered by the UBVRI data, and about 30 times larger 
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than the area covered by the UBVRIJK observations. 
This ensures that the parameter space is properly sam- 
pled and that the errors on the model data are at least 
3, respectively 5, times smaller than the error bars on the 
observed data. 

6. Simulations of the observed data 

The galactic model produces "perfect" stellar objects 
down to a specified magnitude (in fact the model is com- 
puted down to R = 30.5). To build a mock colour cata- 
logue simulating the real data, the magnitudes from the 
galactic model are taken as input and the following steps 
are considered: 

— The error in the magnitude (as computed by 
SExtractor) as a function of magnitude for each band 
(based on the properties of the actual image and real 
data). Additionally, one can scale these errors by a 
factor and add in quadrature an "external error" , sim- 
ulating for example errors in the photometric zero 
point. The "observed" magnitude is then drawn from 
a Gaussian with the estimated error as width and cen- 
tred around the "perfect" magnitude. The scaling fac- 
tor and "external error" used in the present paper are 
given below. 

— The saturation at bright magnitudes. The ADU level 
in the central pixel of an object is simulated, and com- 
pared to the actual value used by SExtractor to flag a 
saturated object. 

— The completeness at the faint end. Artificial star tests 
using the EIS images themselves have not yet been 
performed to estimate the completeness functions. As 
a temporary solution the completeness functions have 
been derived for the i?-band and X-band by comparing 
the EIS data to the deeper observations of respectively 
Wolf et al. (2001) and Saracco et al. (2001) over the 
parts where their observations overlap with the EIS ob- 
servations. For the other bands the shape of the func- 
tion that describes the completeness is kept the same 
(for UBVI taking R as a reference, and for J tak- 
ing K as the reference), and scaling the characteristic 
magnitude where the incompleteness sets in using the 
5cr limiting magnitude of the bands as a reference. For 
the stellar catalogue considered here the exact details 
of the completeness functions are not of too much im- 
portance as the estimated magnitudes where the EIS 
data is estimated to be 50% complete (U = 25.1, B = 
26.6, V = 25.5, R = 25.4, / = 24.3, J = 22.9, K = 
21.2) are at the very faint end of the stellar catalogue. 
In other words, the depth of the stellar catalogue is set 
by the ability to separate point from extended sources 
and the remaining contamination by QSOs and unre- 
solved galaxies, not the depth or incompleteness of the 
observations as such. 

These same observations are used to compare the dif- 
ferences in the magnitudes of the objects in common 
between EIS and the external data. Slices in magni- 



tudes are made, and the average error for the EIS ob- 
servations, the average error for the comparison cat- 
alogue (hence Wolf et al. for R, and Saracco et al. 
for K), and the rms of the difference of the stars in 
common determined. The errors in the external cata- 
logues are always smaller than in the EIS data, as ex- 
pected since these observations are deeper than ours. 
However, both underestimate the error as compared 
to the rms in the difference. Based on these compar- 
isons an external error of 0.01 mag in UBVRI and 
0.03 mag in JK, and a scaling factor of two in all 
passbands are adopted in the simulation of the errors. 
With this description the width of the stellar sequence 
in the colour-colour diagrams is described extremely 
well for all magnitudes and passbands (see Fig. ^). 

The simulation also takes into account the S/N limits 
imposed on the single passband and colour catalogues. 
Furthermore, objects that are saturated in one or more 
bands are removed from the mock catalogue, as masks 
are placed around them in the real data and only objects 
outside all masks are being considered in the final colour 
catalogue. 

Finally, two types of mock stellar catalogues are pro- 
duced. One considering all simulated point sources that 
are potentially "observable", and the second catalogue 
that also takes into account the way that point sources are 
being selected from the real colour catalogue by checking if 
the object is brighter than the morphological classification 
limit used in the production of the observed data set in 
at least one passband (see Tab. |l|). This is a conservative 
estimate, as for the real data the value of the SExtractor 
parameter CLASS_STAR also plays a role. However this fea- 
ture in fact allows us to check if the CLASS_STAR values 
chosen are not too low (resulting in too many galaxies and 
leading to high stellar number counts at the faint end) or 
too conservative (leading to number counts that are below 
the predictions). Details on the simulation of the photo- 
metric errors, the completeness model and the simulated 
datasets will be given elsewhere (Groenewegen et al. 2002, 
in preparation). 

7. Results and analysis 

In this section the observed number counts, colour-colour 
and colour-magnitude diagrams, and colour-distribution 
histograms are presented and compared to the simula- 
tions. The effectiveness of using the x 2 -technique is dis- 
cussed. Finally, catalogues of stellar sources, QSOs, and 
unresolved galaxies are presented. Due to space con- 
straints only representative figures will be shown and the 
full set of figures is only available through the on-line edi- 
tion of A&A. Note that most of the figures need to be 
viewed or printed in colour. 

7.1. Number counts 
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7.1.1. General description 

Figure || shows the number counts for every passband 
from the five passband catalogue for UBVRI and from 
the seven passband catalogue for JK. In the lower right 
panel the i?-band number counts are shown based on the 
seven passband catalogue. The agreement with that from 
the five passband catalogue is excellent, and the same ap- 
plies when comparing the number counts in the other op- 
tical bands as derived from the five and seven passband 
catalogues. 

In the left-hand panels the predictions by the galactic 
model are shown, namely the contributions by the disk 
(dashed and solid thin black histograms), the halo (dashed 
and solid blue histogram) and their sum (dashed and solid 
thick black histogram). Shown are the predictions from the 
"perfect" model (the blue and black dashed histograms), 
and the predictions taking into account the error model, 
saturation, completeness and the MAG_STAR_LIM limits 
that are used to create the colour catalogue (the blue and 
black full histograms). The left-hand panels illustrate the 
relative contributions of disk and halo at various magni- 
tudes. In this particular field the contribution of the disk 
exceeds that of the halo for almost all observed magni- 
tudes. The model predicts, however, that obtaining deeper 
images and assuming one would be able to do the (P/E) 
classification at the level of J, K ss 23 — 24 the number 
counts should be dominated by the halo. In the optical 
the halo dominates the disk only at much fainter mag- 
nitudes (I 25, V » 28). Deeper infrared data could 
therefore better constrain the IMF and other parameters 
describing the halo than the present data. 

In the right-hand panels, the model predictions (again, 
the full thick black histogram) are compared to various 
datasets. The blue histogram with the error bars (based 
on >/N statistics) represents the observed stellar number 
counts from the colour catalogue based on the (P /E) clas- 
sification and subsequent application of the x 2 -techniquc 
described in Sect. f|. The pink histogram represents the 
number counts using the (P /E) classification for that band 
only. As, by construction, there are no stars fainter than 
MAG_STAR_LIM in the single passband catalogue, but this 
value usually does not coincide with the borders of the 
binning, the counts in the last bin may drop artificially. 
Therefore, the counts in the last bin of the pink histogram 
are corrected for this effect by multiplying the counts in 
the last bin by the bin width divided by the actual width 
over which there is data. 

7.1.2. Comparing counts from single passband and 
colour catalogues 

Up to MAG_STAR_LIM, indicated by the vertical dashed 
line, the counts from the single passband and the colour 
catalogue are in very good agreement. In fact, naively, one 
might even have expected them to be identical. They are 
not for several reasons. 
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Fig. 2. Number counts in UBVRI from the five passband 
catalogue and in JK from the seven passband catalogue. 
This sample is for the i?-band, and the full figure is only 
available in the on-line edition. In the lower right panel the 
i?-band counts based on the seven passband catalogue are 
presented. They are in good agreement with those based 
on the five passband catalogue, as is the case for U BVI. 
Left panel: shown are the disk component (dashed and 
solid thin black histograms) and halo component (dashed 
and solid blue histograms), and their sum (dashed and 
solid thick black histograms) predicted by the galactic 
model. Distinguished are the predicted model counts tak- 
ing into account the error model, saturation, completeness 
and the MAG_STAR_LIM in the individual passbands (the 
blue and black full histograms), as well as the "perfect" 
model (the blue and black dashed histograms). The mor- 
phological classification limits for the bands are indicated 
by the vertical dashed line. Right panel: comparison of 
the data (blue coloured histogram with errors bars) with 
the model. Note the different magnitude scale w.r.t. the 
left-hand panel. Shown are the predicted counts taking 
into account the error model, saturation and complete- 
ness (thick dotted black line), as well as the counts that 
additionally take into account the MAG_STAR_lim limits in 
the individual passbands (thick solid black line; the same 
curve as in the left-hand panel). The pink histogram rep- 
resents the observed number counts using the point source 
classification for that band only. The number counts of ob- 
jects assigned "stellar" template spectra that are ranked 1 
and < 3 in the full five and full seven passband catalogue 
respectively, based on the x 2 -technique only, are indicated 
by the dashed blue histogram. The QSO numbers counts 
(among the objects initially classified as point sources) are 
indicated by the red histogram. The morphological classi- 
fication limits for the bands are indicated by the vertical 
dashed line. 

Firstly, in particular in U, at the bright end, the counts 
from the single passband are above the counts from the 
colour catalogue. This is due to the fact that only objects 
which are outside all masks created at the single passband 
level are considered in the colour catalogue. Since the U- 
band is shallow, the single passband catalogue contains 
objects that are not saturated in U, but are saturated 
(hence will have a mask placed around them) in one of 
the other bands, and will not make it into the colour cat- 
alogue. 
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Secondly, the "good" area over which the sources are 
considered is slightly larger for the single passband cata- 
logues compared to the colour catalogue, as the different 
single passband images do not exactly have the same field 
centre and the colour catalogue is created over the area 
that is in common to all passbands. At bright and interme- 
diate magnitudes the counts match to within the Poisson 
errors, as expected (except for the [/-band for reasons just 
discussed). 

Thirdly, when approaching the morphological classifi- 
cation limit the number counts based on the colour cat- 
alogue tend to fall below the counts based on the single 
passband. This is due to the fact that galaxy interlop- 
ers and QSOs have been removed from the sample us- 
ing the x 2 -t ecn nique applied to the objects detected in 
three or more passbands. As illustration, the QSO number 
counts for the objects originally classified as point sources 
is shown by the red histogram. 

Beyond the morphological classification limit, the 
number counts based on the single passband level drop to 
zero (by construction), while the number counts based on 
the colour catalogue extend to significantly fainter mag- 
nitudes. This illustrates the usefulness of using all avail- 
able colours to construct a point source catalogue. This is 
especially true in U, where the counts reach almost two 
magnitudes fainter than the single passband counts. This 
is because the U data is shallow compared to those in the 
other bands and most of the information on the (P/E) 
classification comes from redder passbands. The gain in 
magnitude is negligible in R which, according to Tab. §, 
in most cases defines the (P/E) classification. 

7.1.3. Comparison with the model and literature 

The counts predicted by the model and the observations 
from the colour catalogue are in good agreement, consid- 
ering that no fine tuning or fitting was performed. The ob- 
served data follow the thick black lines within the errors 
indicating that the choice of CLASS_STAR is appropriate 
(or at least consistent with that predicted by the model). 
The data fall below the dashed thick lines, which would 
be the predicted counts if the (P/E) classification would 
be reliable for all magnitudes. Based on the flattening of 
the observed counts and guided by the results from the 
model the 90% completeness limits are estimated to be 
approximately U = 23.8, B = 24.0, V = 23.5, R = 23.0, 
/ = 21.0, J = 20.5, K = 19.0. 

The deepest stellar counts available in the literature 
are from the hst, but they cover a small area, e.g. Santiago 
et al. (1996) combine 17 WFPC-2 fields for a total area 
of 0.027 deg 2 , 90% complete to I 61i of 22 to 24 depending 
on the field. At the other end of the spectrum there is 
the recent analysis of Chen et al. (2001) using 279 deg 2 of 
SDSS data complete to a depth of g' = 21.0 (corresponding 
to approximately V — 20.4 for a typical colour of g' — r' 
of 1, or three magnitudes brighter than ours), or Ojha 
(2001) that analyse seven fields of 2MASS data covering 



67 deg 2 down to K s = 15. Intermediate large area surveys 
are e.g. the "Deep Multicolor Survey" (DMS; Osmer et al. 
1998 and references therein) covering six fields of about 
0.14 deg 2 each complete to about B = 23.5 (comparable 
to ours), or Reid & Majewski (1993) covering 0.3 deg 2 
complete to about V = 21.5 (brighter by two magnitudes 
compared to ours). 

What will make the DPS survey unique when finished 
is that it will cover a larger area, to a depth that is com- 
parable or better, than existing intermediate large area 
stellar surveys and that the data will be in five or seven 
passbands. 

7.1.4. The slope of the IMF 

In this Section the influence of the shape of the IMF on the 
number counts is illustrated. Figure |^ shows the number 
counts in R and J when the power law form in Chabrier 
(2001) is adopted instead of his lognormal law. The num- 
ber counts at R ~ 22 and J ~ 19 — 20 are larger than 
in the default case by about 20%. In the other bands the 
effect is less striking but also present. This illustrates the 
potential of the present data, especially when the DPS sur- 
vey is completed and data of similar depth on fields with 
different galactic coordinates are available, to constrain 
e.g. the IMF. The counts in this field are better fitted 
with the lognormal law proposed by Chabrier (2001) than 
the power law proposed by him. The stars in the relevant 
magnitude interval have typical masses of about 0.1 - 0.2 
M Q . However, the fits are certainly not perfect, and this 
possibly hints to a more complex, or different, shape of 
the IMF than a power law or lognormal law. 

7.2. Colour-colour diagrams 
7.2.1. General description 

Figure [| shows selected colour-colour diagrams split into 
magnitude bins for the five and seven passband catalogues. 
In the left-hand panel the predictions by the model are 
shown. Objects from the halo (open circles) and the disk 
(filled circles), with stars with logg > 7 (representing 
WDs) and stars with masses <0.1 Mq (representing low 
mass stars) additionally marked by a plus and cross sign, 
respectively. The right-hand panel shows the stellar ob- 
jects (filled and open circles), while the objects identified 
as QSOs by the x 2 -technique are shown as crosses. 

The overall agreement is very good in all cases and 
for all magnitude bins. In particular the width of the 
sequences is well matched to the data confirming a- 
posteriori the correctness of scaling the simulated errors 
by SExtractor as was derived by the comparison to deeper 
external data. 

The largest discrepancies are: (a) in the (V — R) — 
(R — I) diagram, where for increasing R — I the model 
predictions saturate at a level of V — R of ~ 1, while the 
observed data points continue to become redder; (b) in 
the (i? — J) — (J — J) diagram, where the slope of the data 
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Number counts in RJ with the power law form of the IMF from Chabrier (2001). The meaning of the curves 
Fig. ||. The counts at R ~ 22, J ~ 19.5 are about 20% larger than for the standard case with the lognormal 



and the model differ for I - J>0.8; (c) in the (V - J) - 
(J — K) diagram, where for V — J > 3 the data seem to 
remain at a constant colour, while there is a slope in the 
model predictions. Possible reasons for the discrepancies 



are discussed in Sect. 7.4.2 



7.2.2. Outliers 

Obviously there are also outliers in the colour-colour di- 
agrams on the data side. These can be real (for example 
due to variability) or hint to problems in the catalogue 
production. This issue is extensively discussed in H2002. 
Most important for the present paper is to remind that 
these outliers do not influence the assignment of the spec- 
tral type as this is based on all available magnitudes and 
errors. 

To guide the eye, outliers are marked by the open cir- 
cles. Outliers are defined in the present paper as follows^]. 
An imaginary box in colour space is drawn around all of 
the model stars in the left-hand panel, for each magni- 
tude bin. An object is defined an outlier if an observed 
data point on the right is not within any of these boxes. 
The width of the box was set after some experimentation 
to be 15-30% of the maximum of the range in both colours 
axis in a given magnitude bin. 

This concept is found useful to estimate the possible 
level of remaining contamination by QSOs and unresolved 
galaxies that have not been identified as such by the \ 2 ~ 
technique. Although the outliers are marked in all colour- 
colour diagrams, the diagrams where QSOs are most easily 
identified based on their difference in colours with respect 
to normal stars are (U — B) — (B — V) in five passbands 
and (U-B)-(B- V), (B-V)-(V- K), (V — J) — 
(J — K) in seven passbands (H2002). Note, however, that 
the colour-colour diagrams indicate that QSOs and stars 
overlap in colour-colour space, and therefore no selection 
based purely on colour-colour diagrams can be made (see 
also H2002). 

The number of outliers in the (U—B) — (B — V) colour- 
colour diagram in five passbands is 36, 28 of which are in 



the region occupied by QSOs, while in the (U — B) — (B — 
V), (B-V)-(V-K), (V-J)-(J-K) colour-colour dia- 
grams in seven passbands there are 63 outliers, 26 of which 
are in the region occupied by QSOs. However since stars 
do overlap with QSOs it can not be assumed automati- 
cally that all outliers are QSOs. The upper limits to the 
remaining contamination by QSOs are therefore <6.3% in 
the seven and <2.3% in the five passband catalogue. In 



In H2002, outliers were defined by considering only the data 
itself, using dissimilarity measures. See H2002 for details. 



Sect. 7.5 an independent estimate of this contamination is 
given. 

7.3. Colour-magnitude diagrams 

Figure|| shows selected colour-magnitude diagrams. In the 
left-hand panel the model predictions are shown with sym- 
bols as introduced in Fig. |]. The right-hand panel shows 
the data. The match is good, in particular the Galactic 
halo (with bluer colours) and the Galactic disk (with red- 
der colours) are readily distinguished. 

One can also notice the particular position occupied 
by the model very low-mass stars (crosses), that are the 
reddest and faintest objects in most plots. White dwarfs 
(plus signs), although found over a relatively wide range 
of colours and magnitudes, are the only objects to occupy 
a vertical strip to the blue of the bulk of halo stars: in 
fact, in diagrams like V vs. V — R and I vs. V — I, the 
sequence of hot WDs (with V - R < 0.3, V — I < 0.4) 
starts at V > 20 and I > 20, and this matches quite well 
the sequence of hot objects seen in the observational plots. 

7.4. Colour distribution 
7.4.1. General description 

Although the colour-colour and colour-magnitude dia- 
grams provide useful and complementary information, the 
most stringent test for the comparison between data and 
model is provided by the distributions over colour at dif- 
ferent magnitude slices, as provided in Figure ^. 

The overall agreement is again very good. The shape 
of the colour distributions is correctly recovered in most 
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Fig. 4. Colour-colour diagram for different cuts in mag- 
nitude. This (V — R) — (R — I) diagram based on the five 
passband catalogue is a typical example and the full figure 
is only available through the on-line edition. Left panel: 
Model predictions. Halo (open circles) and disk (filled cir- 
cles), with stars with \ogg > 7 (representing WDs) and 
stars with masses <0.1 M© (representing low mass stars) 
additionally marked by a plus and cross sign, respectively. 
Since the full simulation is run over n— times the observed 
area, only — of all model data points are plotted. Right 
panel: the data. Open circles represent outliers, in the 
sense that they are not within a certain magnitude dif- 
ference of any of the model data points in the left-hand 
panel (see text). Red crosses are objects identified as QSOs 
by the x 2 -technique. 



cases; the absolute number of objects in the faintest bin, 
however, is not always exact (e.g. U — B and R— I). 

The most obvious discrepancies are: (a) the excess of 
objects with -1 < U — B < -0.5 and 22 < U < 24 which 
could be QSOs not identified by the x 2 -technique, or WDs 
not predicted by the model; (b) the fact that for disk stars 
the predicted V — R colours are bluer than the observed 
ones; (c) similarly, that for disk stars the predicted R — I 
colours are bluer than the observed ones; (d) the excess 
of observed sources with J — K > 1.2 and 18 < K < 19 




0.5 1 1.5 2 0.5 1 1.5 2 



V-R V-R 

Fig. 5. Colour-magnitude diagram. This (V, V — R) di- 
agram based on the five passband catalogue is a typical 
example and the full figure is only available through the 
on-line edition. Left panel: model results. The symbols 
have the same meaning as in Fig. |J. Since the full simula- 
tion is run over n— times the observed area, only — of all 
model data points are plotted. The right panel shows the 
data. 



which again could be QSOs, or galaxies, not identified by 
the x 2 -techniquc. 

The shifts in colour seen in V — R and R — I are of 
the order 0.2 mag, but the conclusion is that all the other 
colours are correctly predicted at a level <0.1 mag. 

7.4.2. Further investigation of the model data 

To make better use of the simulations a numerical code 
has been developed that allows to search on any combi- 
nations of the parameters provided by the Galactic Model 
convolved with the error model, completeness, etc, as de- 
scribed earlier. The parameters provided are the colours 
with their errors, the Galactic component (i.e. disk, halo, 
bulge), age, [Fe/H], initial mass, luminosity, effective tem- 
perature, log g, distance modulus, visual extinction and 
apparent bolometric magnitude. Constraints on the dif- 
ferences between two columns are allowed for as well. As 
an illustrative example the problem of the discrepancy 
between model and data in the (V — R) — (R — I) colour- 
colour diagram is addressed. Figure shows the distribu- 
tion over various model parameters of all simulated stars 
that are formed in the disc and fulfil 1.4 < R — I < 3.5 
and 23 < R < 26. From this analysis it is clear that 
mainly low mass stars (< 0.2 M©) are contributing to this 
colour range indicating that likely the theoretical colours 
for these stars are off by approximately 0.2 mag. 

It should be pointed out that, among the possible in- 
adequacies of the model, systematic shifts in the model 
effective temperatures or in the adopted stellar metallici- 
ties would affect all colours in more or less the same way, 
which does not seem to be the case here. In fact, just a 
few of the colours are deviant, and this might be related 
to the tables of bolometric corrections which are adopted 
in the model. For instance, either these bolometric cor- 
rections were derived from imperfect response curves for 
some of the EIS filters, or the synthetic stellar spectra in 
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Fig. 6. Colour distributions for different cuts in magni- 
tude. The full figure is only available through the on-line 
edition. Shown are the model predictions (the halo is the 
blue lined histogram; the disk the thin lined histogram, 
and their sum is the thick lined histogram) and the data 
(the pink histogram with the error bars). 

use (see Girardi et al. 2002 for details) present inadequa- 
cies in some of their wavelength regions (caused by e.g. 
incomplete line opacity tables). This matter will be fur- 
ther investigated, but whatever is their cause, these colour 
problems are small enough to not affect any of our con- 
clusions. 

Alternatively, these discrepancies may result from the 
data rather than the model. It is important to point out 
that most of the discrepancies occur at the near-infrared 
passbands (UK) where the data may be affected by fring- 
ing as well as other problems related to the reduction of 
jittered infrared data. It will be of interest to see if these 
problems persist as the techniques for reducing WFI and 
SOFI data are improved. 

7.5. The effectiveness of the x 2 -technique 

In the selection of stellar sources, the values of the 
CLASS_STAR of the different passbands are combined as 
described previously. Then the x 2 -technique is used to fil- 
ter out unresolved galaxies and QSOs. In this section it is 
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Fig. 7. Distributions over various parameters of all model 
stars which are formed in the disc and have 1.4 < R — I < 
3.5 and 23 < R < 26. 



investigated how reliably the % 2 -technique works, and if 
it could be used as a stand-alone tool. 

The dashed blue histogram in Fig. || in fact shows the 
number counts of objects that have been assigned "stellar" 
template spectra and that are ranked 1 and < 3 when 
applying the x 2 -technique to the full, five, respectively 
seven, passband colour catalogue, not just the point source 
catalogue^. 



A technical remark: The difference of selecting objects 
ranked 1 for the five passband catalogue and ranked <3 (i.e. 
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Even at relatively bright magnitudes these counts are 
significantly larger than predicted by the model, in partic- 
ular in UBVR. This is very likely due to misclassifications, 
and cross-talk, between galaxies and stars. Since the num- 
ber density of galaxies is so much larger than that for stars 
(already a factor of 10-15 at mag_star_lim), even a 10% 
misclassification of star to galaxies and vice versa leads 
to a very small decrease of galaxy number counts while it 
leads to a doubling of the stellar number counts. This is 
illustrated further below, and shows that the x 2 -technique 
can not, and should not, be used blindly. 

Fig. H shows the distribution of CLASS_STAR in the 
different bands for the objects that are classified point 
sources in the colour catalogue (red circles). This gives an 
illustration of the usefulness of using all available colours 
to classify objects. In UV many objects classified as point 
sources using all available colours are outside the region 
that is initially set to contain point sources in these bands 
individually. Table |^ contains information on which band 
actually decided whether an object is classified as a point 
source. This is B, R and / for the five passband catalogue, 
and B, R and J for the seven passband catalogue. The dots 
in these figures show the distribution of CLASS_STAR for 
the objects assigned "stellar" template spectra that are 
ranked 1 and < 3 in the full five and full seven passband 
catalogue respectively, when applying the x 2 ~technique. 
It can be clearly seen that even bright objects with small 
CLASS_STAR values can be assigned a stellar spectral type, 
even when these objects are verified to be resolved galax- 
ies. 

This illustrates the strength and weakness of using 
MAG_STAR_lim/classj3TAR_lim or the x 2 -technique. The 
former method selects point like objects without having 
the possibility to distinguish between stars and QSOs or 
unresolved galaxies, while the latter method uses astro- 
physical information (the SEDs of objects in the Universe) 
to assign types but it does not consider whether an object 
is resolved or not. 

To investigate the level of misclassification of true 
galaxies as stars by the x 2 -technique, Fig. ^ shows the 
ratio of the number of objects assigned "stellar" spectra 
to all objects, for values of class_star < 0.1 (i.e. true 
galaxies). It is noticed that the fraction of misclassifica- 
tion is essentially independent of magnitude and ranges 
between 13 and 19% in the five passband catalogue (with 
an average of 16.3%), and between 7 and 17% for the 
seven passband catalogue (with an average of 12.8%). The 
smaller level of contamination indicates the effect of hav- 

all objects) for the seven passband catalogue is related to the 
information provided by the extra filters. More precisely, for 
objects that have an infrared continuum described by a power 
law, such as low- to intermediate redshift QSOs or hot WDs, 
the J and K magnitudes do not provide any additional infor- 
mation about the SED but will impose additional constraints, 
by increasing the number of the degrees-of-freedom in the fit- 
ting. A possible answer to the problem of redundant informa- 
tion could be a Principle Component Analysis. This will be 
investigated in the future. 



ing more colours available in the fitting to better constrain 
the assignment of the best fitting spectrum. 

If it is assumed that the misclassification from true 
stars to galaxies or QSOs is similar, one can estimate from 
the numbers in Table | that 16.3% out of the total of 173 
QSO+unresolved galaxies, or 28 objects may be wrongly 
assigned such a type, or that 28 stars are incorrectly re- 
moved from the stelar catalogue. This constitutes 2.4% of 
the sample of stars. The same number is found for the 
seven passband catalogue. 
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Fig. 8. CLASS^STAR values of the objects that are classi- 
fied as point sources in the five and seven passband colour 
catalogues (filled red circles). This panel shows the i?- 
band based on the five passband catalogue as an example, 
and the full figure is only available through the on-line 
edition. Black dots represent the CLASS_STAR values of 
the objects that are classified as "stellar" and that are 
ranked 1 and < 3 in the full five and full seven passband 
catalogue respectively, when applying the x 2 ~technique. 
mag_STAR_lim and CLASSJ3TAR limits are indicated. 



7.6. Spectral types and data tables 

Tables | and | give the first entries of the stellar catalogue 
in respectively five and seven bands. The tables give the 
following information: in column (1) the EIS identification 
number; in columns (2) and (3) the J2000 coordinates; in 
columns (4) - (13) for the five passband and columns (4) - 
(17) for the seven passband the magnitude (Vega system) 
and the SExtractor error (upper limits are indicated by a 
— 1; in column (14) (respectively column (18) in the seven 
passband) the assigned spectral type, with the following 
meaning: 

— Prefix "MS" . Templates for normal O to M stars from 
the Pickles (1998) library. Stars that are detected in 
fewer than three bands and therefore have not been 
fitted by the x 2 -technique are listed as "dummy" here. 
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Fig. 9. Fraction of objects that are classified as "stellar" 
and that are ranked 1 and < 3 in the full five and full 
seven passband catalogue respectively, when applying the 
X 2 -technique and that have a CLASS_STAR less than 0.1. 
MAG_STAR_LIM limits are indicated. This panel shows the 
i?-band based on the five passband catalogue as an exam- 
ple and the full hgure is only available through the on-line 
edition. 

— Prefix "WD". Templates for WDs from the models 
provided by D. Koster (the numbers indicate effective 
temperature and logg), Ibata et al. (2000; the observed 
spectrum of F351-50, F821-07) or Oppenheimer et al. 
(2001; the observed spectrum of WD 0346+246). 

— Prefix "LMS" . Templates for Low Mass Stars from 
Chabrier et al. (2000). The prefix is followed by a num- 
ber indicating T c ff/100 and the type of model used 
(NeGen = Next Generation, AMESd = DUSTY, or 
AMESc = CONDENSED, as defined in Chabrier et 
al. 2000). 

Tables |[ @ and ||give the first entries of the likely QSOs 
and unresolved galaxies in the five and seven passbands 
as identified by the procedure in Sect. ^. The information 
listed is the same as in the previous tables, except that in 
the last column the photometric redshift is listed for the 
QSOs, and the type of galaxy (following Coleman et al. 
1980) and photometric redshift for the unresolved galaxies 
as determined by the spectral template fitting method. 

The complete tables can be re- 
trieved from the CDS or from the URL 



" |http://www. eso.org/science/eis/eis_pub/cis_pub. htm] 



8. Summary and conclusions 

This paper describes the procedures adopted in the con- 
struction of multi-colour stellar catalogues, suitable for 
statistical studies, extracted from multi-band imaging 
data. The procedure involves several steps among which: 
a magnitude-dependent scheme to morphologically clas- 
sify sources identified using SExtractor; a procedure to 



combine the information available in catalogues extracted 
from images taken in different passbands; and the use of a 
% 2 -technique to assign spectral types to the sources using 
the multi-band information. This allows to minimise the 
contamination of the stellar catalogue by QSOs and unre- 
solved galaxies, and to assign spectral types for the stars, 
which are included in the stellar catalogues presented in 
this paper. 

The methodology outlined here is applied to the first 
data set released for the DPS, the CDF-S. The 90% com- 
pleteness limits are estimated to be approximately U — 
23.8, B = 24.0, V = 23.5, R = 23.0, / = 21.0, J = 20.5, K 
= 19.0. To assess the quality of the catalogues produced, 
number-counts, colour-colour and colour-magnitude dia- 
grams, and colour-distributions derived from the data are 
compared to those obtained from simulated catalogues. 
Mock catalogues are created using a Galactic model based 
on population synthesis (described by parameters set by 
independent data), an error model describing the expected 
photometric errors, saturation and completeness as de- 
rived from the data. Even though no attempt is made to 
fine tune the Galactic model parameters the agreement 
between real and simulated data is remarkable, serving as 
a good indicator of the reliability of the catalogues pro- 
duced. The comparison also suggests that the depth of 
the data is suitable to constrain the IMF and/or SFR of 
low-mass stars. 

This paper represents a first attempt to define proce- 
dures to produce well-defined, deep stellar catalogues with 
minimal contamination by QSOs and unresolved galaxies. 
The results are encouraging and demonstrate the valu- 
able contribution that the homogeneous multi-band opti- 
cal/infrared data set from DPS, probing different direc- 
tions of the Galaxy, will make when completed. 
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Table 3. First 15 entries of the stellar catalogue in five passbands. 
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Table 4. First 15 entries of the stellar catalogue in seven passbands. 
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Table 5. First 15 entries of the likely QSOs in five passbands. 
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Name 




a[ JzOOO J 


o[ J2000 j 


l, ! 


0"T; 


B 


0"B 


V 




R 




I 


m 


.7 




R 


0"I< 


2 phot 


EISJ033145.22- 


275435.8 


03:31:45.22 


-27:54:35.8 


19.32 


0.00 


20 


04 





00 


20.01 


0.00 


19.81 


0.00 


19.37 


0.01 


19.72 


0.08 


18.09 


0.04 


0.50 


EISJ033145.50- 


275446.4 


03:31:45.50 


-27:54:46.4 


22.31 


0.02 


22 


88 





01 


22.63 


0.03 


21.85 


0.02 


21.17 


0.02 


20.67 


0.13 


19.05 


0.07 


0.20 


EISJ033147.98- 


275045.5 


03:31:47.98 


-27:50:45.5 


23.18 


0.05 


23 


63 





02 


23.78 


0.06 


23.42 


0.05 


22.85 


0.07 


21.76 


0.1 1 


20.25 


0.13 


2.90 


EISJ033148.26- 


274330.2 


03:31:48.26 


-27:43:30.2 


22.40 


0.03 


23 


25 





02 


22.71 


0.03 


22.14 


0.02 


22.51 


0.05 


22.14 


0.13 


21.19 


0.30 


0.10 


EISJ033149.43- 


274634.4 


03:31:49.43 


-27:46:34.4 


22.22 


0.02 


22 


84 





01 


22.44 


0.03 


22.21 


0.02 


21.92 


0.04 


20.81 


0.05 


19.17 


0.06 


0.65 


EISJ033150.97- 


274116.2 


03:31:50.97 


-27:41:16.2 


22.23 


0.02 


22 


86 





01 


22.66 


0.03 


22.02 


0.02 


21.20 


0.02 


20.04 


0.03 


18.67 


0.04 


1.10 


EIS J033151.20- 


274852.3 


03:31:51.20 


-27:48:52.3 


22.58 


0.0 1 


23 


56 





05 


23.10 


0.04 


22.77 


0.03 


21.76 


0.03 


20.91 


0.06 


19.66 


0.08 


0.20 


EIS J033153.52- 


275541.3 


03:31:53.52 


-27:55:41.3 


23.28 


0.05 


23 


56 





02 


23.42 


0.05 


22.90 


0.03 


23.01 


0.07 


21.91 


0.17 


20.04 


0.1 1 


2.94 


EISJ033155.63- 


275402.4 


03:31:55.63 


-27:54:02.4 


22.67 


0.04 


23 


29 





02 


22.93 


0.03 


22.54 


0.03 


22.14 


0.04 


20.76 


0.1 1 


19.74 


0.09 


0.00 


EIS J033157.58- 


275050.4 


03:31:57.58 


-27:50:50.4 


22.84 


0.0 1 


23 


53 





02 


23.30 


0.04 


22.51 


0.02 


21.96 


0.04 


21.32 


0.06 


19.91 


0.09 


0.1 1 


EIS J033158.29- 


275426.5 


03:31:58.29 


-27:54:26.5 


22.16 


0.02 


22 


93 





01 


22.83 


0.03 


22.56 


0.02 


21.89 


0.03 


21.28 


0. 13 


19.87 


0.10 


0.72 


EISJ033158.57- 


273912.9 


03:31:58.57 


-27:39:12.9 


22.68 


0.03 


23 


59 





02 


23.49 


0.05 


23.11 


0.04 


22.89 


0.06 


21.67 


0.15 


21.03 


-1.00 


0.50 


EISJ033200.61- 


274108.0 


03:32:00.61 


-27:41:08.0 


22.69 


0.03 


23 


52 





02 


23.38 


0.04 


22.53 


0.02 


21.92 


0.03 


21.40 


0.07 


20.22 


0.15 


0.40 


EISJ033201.32- 


274006.1 


03:32:01.32 


-27:40:06.1 


22.09 


0.02 


22 


89 





01 


22.79 


0.03 


22.50 


0.02 


22.06 


0.04 


21.04 


0.06 


20.13 


0.12 


0.72 


EISJ033202.25- 


274248.8 


03:32:02.25 


-27:42:48.8 


22.91 


0.03 


23 


76 





02 


23.56 


0.05 


23.10 


0.03 


22.08 


0.03 


20.97 


0.06 


20.49 


0.15 


0.20 



Table 7. First 15 entries of the likely unresolved galaxies in five passbands. 



Name 






a 


(J2000) 


<S(./2000) 




U 


eru 


B 




V 




R 


era 


/ 




Type 


-2phot 


EIS J033114 


30- 


-274707 


5 


03 


31 


14 


30 


-27 


47 


07 


5 


22 


98 


0.06 


23 


68 





03 


23 


06 





05 


22 


39 





03 


21 


88 





04 


Im 


0.4 


EIS J033115 


17- 


-274453 


4 


03 


31 


15 


17 


-27 


44 


53 


4 


23 


43 


0.06 


23 


70 





03 


23 


14 





04 


22 


16 





02 


21 


33 





02 


Sbc 


0.3 


EIS J033117 


08- 


-274617 


9 


03 


31 


17 


08 


-27 


46 


17 


9 


24 


93 


0.16 


24 


05 





03 


22 


94 





04 


21 


53 





01 


20 


29 





01 


E 


0.4 


EIS J033121 


10- 


-274043 


2 


03 


31 


21 


10 


-27 


40 


43 


2 


22 


49 


0.03 


23 


36 





02 


23 


01 





04 


22 


38 





02 


21 


33 





03 


Sbc 


1.6 


EIS J033121 


74- 


-275943 


1 


03 


31 


21 


74 


-27 


59 


43 


1 


22 


16 


0.02 


22 


93 





01 


22 


68 





03 


21 


90 





02 


21 


27 





02 


SB3 


0.5 


EIS J033122 


79- 


-275932 


9 


03 


31 


22 


79 


-27 


59 


32 


9 


22 


89 


0.03 


23 


46 





02 


22 


96 





03 


22 


32 





02 


22 


02 





04 


Tm 


0.3 


EIS J033128 


52- 


-275135 


9 


03 


31 


28 


52 


-27 


51 


35 


9 


22 


56 


0.03 


23 


1 1 





01 


22 


65 





03 


21 


67 





01 


20 


84 





02 


Scd 


0.5 


EIS J033131 


14- 


-273343 


5 


03 


31 


31 


14 


-27 


33 


43 


5 


23 


44 


0.06 


23 


64 





02 


22 


79 





03 


22 


07 





02 


21 


47 





03 


Sbc 


0.3 


EIS J033131 


14- 


-275846 


8 


03 


31 


31 


14 


-27 


58 


46 


8 


19 


24 


0.00 


20 


00 





00 


19 


82 





00 


19 


41 





00 


18 


14 





00 


E 


2.7 


EIS J033131 


27- 


-273429 


6 


03 


31 


31 


27 


-27 


34 


29 


6 


23 


13 


0.04 


23 


46 





02 


23 


48 





05 


22 


73 





03 


21 


64 





03 


E 


2.8 


EIS J033132 


23- 


-275600 


5 


03 


31 


32 


23 


-27 


56 


00 


5 


22 


36 


0.02 


23 


31 





02 


23 


29 





04 


22 


81 





03 


21 


48 





03 


Scd 


1.4 


EIS J033134 


56- 


-275647 


2 


03 


31 


34 


56 


-27 


56 


47 


2 


22 


70 


0.03 


23 


46 





02 


23 


1 1 





04 


22 


33 





02 


22 


03 





03 


SB2 


0.4 


EIS J033143 


67- 


-274645 


4 


03 


31 


43 


67 


-27 


46 


45 


4 


22 


73 


0.04 


23 


05 





01 


22 


46 





03 


21 


54 





01 


20 


81 





02 


Scd 


0.4 


EIS J033146 


67- 


-274126 


8 


03 


31 


46 


67 


-27 


41 


26 


8 


25 


28 


-1.00 


25 


04 





07 


23 


58 





06 


22 


32 





02 


20 


98 





02 


E 


0.5 


EIS J033149 


55- 


-274319 


7 


03 


31 


49 


55 


-27 


43 


19 


7 


22 


50 


0.03 


23 


04 





01 


22 


40 





02 


21 


69 





01 


21 


21 





02 


Scd 


0.3 



Table 8. First 15 entries of the likely unresolved galaxies in seven passbands. 



Name q(J2000) <5(J2000) U <ru B <r B V <r v R <tr / a\ J <tj K <r K Type z pho i 

EIS J 033148. 76— 274250.3 03:31:48.76 -27:42:50.3 25.25 0.21 25.02 0.06 24.33 0.08 22.57 0.03 20.78 0.02 19.40 0.02 17.93 0.02 E 

EIS J033150. 32— 274241.1 03:31:50.32 -27:42:41.1 25.28 -1.00 25.42 0.07 24.16 0.07 22.47 0.02 21.09 0.02 19.81 0.02 18.37 0.03 E 

EIS J033159. 69-274701.6 03:31:59.69 -27:47:01.6 25.28 -1.00 25.29 0.09 24.98 0.20 23.10 0.06 21.24 0.03 19.67 0.06 17.70 0.03 E 1 

EIS J033216. 21-273930.6 03:32:16.21 -27:39:30.6 19.92 0.00 20.51 0.00 20.20 0.01 19.65 0.00 19.23 0.01 18.70 0.02 17.53 0.02 SB2 

EISJ033217.81-275620.9 03:32:17.81 -27:56:20.9 22.14 0.02 22.90 0.01 22.68 0.03 21.99 0.02 21.10 0.02 20.38 0.04 19.17 0.05 SB3 

EIS J033219. 89-274721.4 03:32:19.89 -27:47:21.4 25.28 -1.00 25.65 0.09 23.95 0.09 22.60 0.03 21.16 0.02 19.98 0.02 18.19 0.02 E 

EIS J033221. 67— 275502.0 03:32:21.67 -27:55:02.0 25.28 -1.00 25.00 0.05 23.93 0.07 22.52 0.02 21.29 0.02 20.27 0.04 18.62 0.03 E 

EIS J033222. 14— 275101.9 03:32:22.14 -27:51:01.9 25.28 -1.00 25.58 0.09 24.16 0.08 22.92 0.03 21.40 0.02 20.10 0.02 18.45 0.02 E 

EIS J033225. 93-274401.6 03:32:25.93 -27:44:01.6 22.29 0.02 22.84 0.01 22.33 0.02 21.86 0.02 21.50 0.02 20.99 0.05 20.02 0.08 SB2 

EIS J033239. 24-274758.7 03:32:39.24 -27:47:58.7 25.28 -1.00 24.54 0.06 23.24 0.05 21.79 0.02 20.29 0.01 19.06 0.01 17.41 0.01 E 

EIS J033239. 48— 275032.0 03:32:39.48 -27:50:32.0 25.65 0.30 24.84 0.05 23.73 0.05 22.16 0.02 20.89 0.02 19.65 0.01 18.00 0.03 E 

EIS J033241. 77-275101.0 03:32:41.77 -27:51:01.0 25.28 -1.00 25.00 0.05 23.89 0.06 22.41 0.02 21.27 0.02 20.11 0.02 18.56 0.04 E 

EIS J033242. 04-275226.2 03:32:42.04 -27:52:26.2 22.66 0.03 23.37 0.02 23.23 0.04 22.32 0.02 21.81 0.03 20.99 0.04 19.91 0.09 SB3 

EIS J033251. 16-274431.9 03:32:51.16 -27:44:31.9 25.28 -1.00 24.53 0.04 23.54 0.06 22.16 0.02 20.94 0.02 19.82 0.02 18.28 0.03 E 

EIS J033251. 67-273937.1 03:32:51.67 -27:39:37.1 23.11 0.04 23.70 0.02 23.45 0.04 23.07 0.03 22.31 0.04 20.66 0.06 19.76 0.09 SB3 1 
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Appendix A: Determining mag_star_lim and 

CLASS_STAR_LIM 



< 
E-i 
m 

in 
GO 
< 




mm 

mm 



20 

R— Magnitude 



2', 




20 

R— Magnitude 

Fig. A.l. CLASS_STAR as a function of i?-band magnitude, 
with the binning in magnitude and CLASS_STAR indicated, 
for the first pass (top), and second pass (bottom, without 
data points for clarity). The hatched bin is the one that 
defines the background of galaxies. 



This appendix gives a brief technical description of 
the way CLASS_STAR_lim and MAG_STAR_LIM are deter- 
mined in the present paper. Refinements of this scheme 
might be considered in the future after completion of the 
DPS survey (12 times the area under consideration here in 
different directions) and analysis of their respective stel- 
lar catalogues, spectroscopic follow-up studies that indi- 
cate the (in)correctness of the assigned spectral types, or 
more extensive numerical simulations including galaxies 
and QSOs. 

One way to analyse a diagram like Fig. |l| is to set 
two levels of significance. The first, ti, setting, for a fixed 
magnitude range, the level of significance of detecting 
stars against the background of galaxies at some value 
of CLASS_STAR, and t<i , setting the level of significance of 



detecting stars at some magnitude, given a fixed range in 

CLASS_STAR. 

This idea is implemented in the following way. The first 
step is to allow for the single passband catalogue to be cut 
at a magnitude error level, and keeping only objects with 
"good" SExtractor flags. 

For the present analysis, only objects with SExtractor 
flag <4 are kept and the default cut of S/N > 2 is imposed. 
The objects retained are divided in bins of magnitude and 
CLASS_STAR. Instead of using a fixed bin size in magnitude 
and CLASSJ3TAR the bins are constructed in such a way 
that they contain an equal number of objects, that is, the 
width of the bin becomes the variable, rather than the 
number of objects in a bin of fixed size. The number of 
bins used is set automatically and is based on the total 
number of objects. 

This is illustrated in Fig. A.l (top panel) for the 
i?-band, where the objects are divided into 14 bins in 
magnitude and 53 bins in CLASS_STAR in the first pass. 
When the density of objects is high the individual bins in 
CLASS_STAR are no longer discernible. 

For each magnitude bin, the distribution over 
CLASS_STAR is analysed in the following way. The bin 
with the largest width in CLASS_STAR is considered to de- 
fine the background of galaxies (indicated by the hatched 
area in Fig. A.l). A correction is made for the number 
of galaxies relative to the total number of objects in that 
bin based on the ratio of the number of objects below 
this bin to the total number of objects. For simplicity 
it is assumed that the background level is constant for 
the bins at higher CLASS_STAR. The mean and rms (based 
on \[N statistics) per unit CLASS_STAR are determined. 
This is a necessary step as all bins have a different width. 
For the bins that have a CLASS_STAR value higher than 
that defining the background bin, the background level 
expected for that bin width is subtracted and divided by 
the rms level expected for that bin width. This gives a 
significance level per bin, Chin- For an external threshold 
value, ii, the lowest bin that has a significance above this 
threshold is determined and, by linear interpolation using 
the next highest bin, the final value for CLASS_STAR_LIM 
for that magnitude bin is calculated. In addition, all the 
bins that have a CLASS_STAR larger than the bin that de- 
fines Chin are combined and a similar significance level for 
the group of stars as a whole (cr groU p) is computed. Once 
the loop over the magnitude bins is completed, and given 
a second external threshold level, t2, that magnitude bin 
is determined where t% > <x group . By linear interpolation 
using the next brightest magnitude bin, the final value for 
MAG_STAR_LIM is computed. Additionally, the magnitude, 
m cu t, is determined below which stars can not be recog- 
nised with any confidence. Then, for a bin in magnitude 
to the left of mag_star_lim (since mag_star_lim signi- 
fies a right-sided cut-off) and containing an equal num- 
ber of objects as before, the binning in CLASS_STAR and 
the analysis to derive the value for CLASS_STAR_lim is re- 
peated. This whole procedure is done twice, the first time 
using all objects that have passed the selection on mag- 
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nitude error and SExtractor flag, and a second pass only 
retaining objects brighter than 77i cut . This allows for a 
better sampling over the magnitude interval where stars 
can be identified. The binning for the second pass and 
the effect of only retaining objects brighter than m cut 
is illustrated in the bottom panel of Fig. A.f . Figure [l] 



shows for the i?-band the distribution over CLASS_STAR 
and CLASS_STAR_lim and MAG_STAR_LIM for the adopted 
choice of t% — 30 and t\ = 80. 

An independent check on the scheme presented, is to 
consider galaxy counts as the complement to the objects 
selected to be point sources, the argument being that 
if one is too liberal in the choice of MAG_STAR_LIM and 
CLASS_STAR_LIM to classify stars, the remaining number 
of "non-stars" could be inconsistent with known galaxy 
counts. In B, R and K the number of observed objects 
in the colour catalogue not classified as stars in the 
0.5 mag bin below the faintest derived MAG_STAR_LIM 
limit (i.e. 23.35 < B < 23.85, 22.25 < R < 22.75, 
18.30 < K < 22.80) are, respectively, 2209, 1935 
and 389. The number density of galaxies quoted in 
the literature is between 6500-11000, 5000-9000 and 
4000-8000 per deg 2 per 0.5 magnitude bin in BRK 



respectively (see the extensive compilation at http://star 



www. d ur.ac.uk / ~nm / pubhtml / counts / counts.html for 



detailed references). The predicted number of galaxies 
for the effective area is 1710-2893, 1315-2367 and 371-742 
in BRK respectively, in agreement with the numbers 
observed, and indicating that the choice of MAG_STAR_LIM 
and CLASS_STAR_LIM to classify stars is not too liberal. 



